Unconventional Gravitational Excitation of a Schwarzschild Black Hole 
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Besides the well-known quasinormal modes, the gravitational spectrum of a Schwarzschild black 
hole also has a continuum part on the negative imaginary frequency axis. The latter is studied 
numerically for quadrupole waves. The results show unexpected striking behavior near the alge- 
braically special frequency Q = — Ai. This reveals a pair of unconventional damped modes very 
near Q, confirmed analytically. 
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Black-hole linearized gravitational waves propagate on 
the spacetime around the event horizon. Because waves 
escape to infinity and into the event horizon, the system 
is dissipative, and described by its nontrivial spectrum 
in the lower-half frequency plane. The latter consists of 
well-studied quasinormal modes (QNMs) [1-4], and a cut 
along the negative imaginary axis (NIA) giving rise to a 
continuum, about which little is hitherto known. 

The Schwarzschild black hole is the simplest compact 
object in relativity, giving its gravitational excitation 
spectrum (part of which should soon be observable [5]) 
a fundamental status akin to the hydrogen problem in 
quantum mechanics. Interest in this (classical) spectrum 
is surging [6], as it seems to offer clues to the quantum 
theory, in particular to a calculation of the Bekenstein 
entropy in loop quantum gravity and to the quantum 
of area. A correspondence between classical frequencies 
and quantum excitations should not be surprising — cf. 
the harmonic oscillator. We shall be content to have 
these developments serve as context and motivation. 

In this article, the continuum spectrum is evaluated 
and characterized. It is found to oscillate in Imw (see 
Fig. 2 below), with a half-period equal to the spacing [7] 
of the string of QNMs parallel and close to the NIA that 
are studied in Ref. [6]. In particular, for the principal case 
of waves propagating from a source near the horizon to a 
distant observer, the spectrum is sharply dominated by a 
dipole contribution, revealing a pair of QNMs lu± on the 
unphysical side of the cut, very close to the algebraically 
special frequency f2 [8]. This represents unconventional 
damped modes. Their effect is more than indirect: as the 
black hole is given some rotation, oj± seem to cross the 
NIA and emerge onto the physical sheet. This follows 
from precise and nontrivial numerical studies, combined 
with analytics at and near f2. 

For a black hole of mass M (c = G = 2M = 1 below) 
and each angular momentum £, the radial functions tp of 
axial gravity waves are governed by a generalized Klein- 
Gordon or so-called Regge- Wheeler equation 



[dl+uj 2 - V(x)]ip(x,(j) = 



(1) 



[£(£+l)r~ 2 — 3r~ 3 ] accounts for the Schwarzschild back- 
ground [9]. We impose outgoing- wave conditions 
(OWCs) ^(x— >±oo, lu) ~ e %ul \ x \ . (For x — > — oo, waves 
thus go into the horizon.) 

Polar waves obey the Zerilli equation [10] — a Klein- 
Gordon equation with V(x) = V(x) + 2d x W(x), with 



W(r) = r + 



3(r-l) 



r 2 [(£-l)(£+2)r + 3] 



(2) 



x = r + ln(r— 1) is the tortoise coordinate and r the 
circumferential radius. The potential V(r) = (1— r _1 ) x 



T = (£+2)l/[6(£-2)\] = in. Below, we take I = 2, the 
most important case, so T = 4. The solutions ip of (1) and 
ip of the Zerilli equation are related by "intertwining" or 
supersymmetry: tp(x,uj) = [d x +W(x)]i/j(x, uj) [1, 11, 12]. 
Thus, also the two continua are closely related, and we 
focus on (1). 

At f2, (1) has an exact solution g(x, — f2) oc 
exp{-J x dyW(y)} [13]; cf. below (4) for g(x,u). The 
nature of f2 has sometimes been controversial, with both 
a QNM [14] and — in apparent contradiction — a total- 
transmission mode [15] having been claimed to occur. 
Moreover, the numerical spectra near il for small Kerr- 
hole rotation a point to intricate behavior for a j [16]. 
The issues were resolved in Ref. [17]: (a) At il, (1) 
and the Zerilli equation cannot be identified, for the su- 
persymmetry transform relating them is singular there. 

(b) Even for each separate sector (axial or polar), the 
OWCs at are highly subtle [cf. below (4)]; only with 
precise definitions can the mode situation be elucidated. 

(c) For o / 0, f! splits into distinct QNM and total- 
transmission-mode multiplets; still, questions remain, to 
which we will return in the Discussion. This paper vin- 
dicates and augments these findings by uncovering the 
pair of nearby unconventional QNMs lo±. Besides be- 
ing of fundamental significance, this further explains the 
problem's potential for numerical artifacts unless special 
precautions are taken. First, some preliminaries must be 
presented. 

The Green's function G(x,y;t) = f(du/2ir) e~ Mt x 
G(x,y;uj) relates i(>(x,t>0) to {ip(y,0),ip(y,0)}. Clos- 
ing the contour in the lower- half cj-plane (Fig. 1, line a) 
splits G into (a) the large semicircle, giving a prompt sig- 
nal propagating directly from y to x and vanishing after 
a finite t [18, 19]; (b) the QNM poles, causing a ringing 
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FIG. 1: Fourier inversion of G; crosses denote poles, (a) Con- 
ventional contour involving a cut on the NIA. (b) Modified 
contour detouring around the unconventional pole near 
(contour and pole on unphysical sheet shown with broken 
lines). For simplicity the mirror contour is not shown. 




FIG. 2: Plot of g( 7 ). 



signal dominating at intermediate t [20, 21]; and (c) our 
main focus, the branch cut on the NIA: 

AG(x, y; -ij) = G+ (x, y; -fy) - G_ (x, y; -ij) , (3) 

where G±(— i-f) — lim e j G(— ij±e) are continuations 
from ±u> > 0. The physical sheet for G+ (G_) thus lies 
to the right (left) of the NIA. The continuum is given by 
AG, and for 7 J, causes the late-i behavior [18, 20]. 
In general, the Green's function G(x,y) — G(y,x) is 

G(x,y;uj) = J^ 1 (uj)f(y,uj)g(x,uj) , y < x , (4) 

where / (<?) solves (1) with the left (right) OWC, and J = 
gf — fg' is their Wronskian. We define f(x— 00, u>) ~ 
1 • e~ lux and g(x— >oo, u>) ~ 1 ■ e' luJX . At a zero of J, / oc g 
satisfies both OWCs and defines a QNM. 

If V has its support in say [— d, d], the OWCs can be 
imposed at ±d, so the wave equation is integrated over 
a finite length; hence, f,g are analytic in u. If how- 
ever V(x^— 00) ~ ^2 k Vke xkx (here A = 1), Born ap- 
proximation shows that / typically has poles (anomalous 
points) at ui n = —inX/2. These are removable by rescal- 
ing f(u) 1 ^ (w—uj n )f(w), leaving G unaffected. How- 
ever, {ffc} could "miraculously" leave some f(u> n ) non- 
singular. For (1), n — 2T is the only miraculous point 
(for any £), which can be studied analytically [17]. 

Cut strength. For x — ► +00, the centrifugal barrier 
does not scatter [18], so consider (7E is Euler's constant) 



fi In r f°° 
V ( x ) ~ 4 ~ -T = / dX (3- 2 7E- In A) A 

x x J Q 



2 e -X X 



(5) 



The superposition of tails e~ Xx spreads the poles — in\/2 
into a cut. On the NIA, g± ~ 1 • e JX obey the 
same wave equation (1). Hence Ag ~ • e lx [cf. (3) 
for A] is the small solution oc g(+if). Since g(— u*) = 
g*{uj), one has ReAg = 0. The cut strength q(-f) = 
Ag(x, — ij)/[ig(x, [17] is independent of V(x) at 

any finite x: if, say, V\(x>L) = V2(x>L), then qi = qi- 
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FIG. 3: Plot of -iAG(x,y; -i-y) for x = 0.2, y = 0.1. 



Thus, q efficiently characterizes the cut in g, causing the 
one in G; cf. (7). Ref. [20] has an expression for q, but it 
is hard to evaluate. 

We integrate g(x, +i"f) from large x, but for g(x, —i"f) 
this would be unstable [22]. Instead, we find g(x, —ij±e) 
by Miller's algorithm of downward recursion on Leaver's 
series [23-25], taking e | in the difference for Ag 
(Miller's algorithm fails at e = 0). Fig. 2 shows the en- 
suing q(j), consistent with q'(0) — 2ir [18]. This and 
subsequent figures constitute an independent look at the 
Schwarzschild spectrum down in the complex aj-plane, 
now thought to have deep physical significance [6] . 

The zeros at 7 = 0.75, 1.35, 1.90, 2.44, 2.96, 3.48, 
4.00, 4.51, 5.03, 5.54, . . . suggest that the spacing ap- 
proaches 5. Indeed 5(7) ~ 4cos(27T7) + 0(7~ 1//2 ) by 
a WKB analysis [26]. The zero at T and moreover 
q'(T) w -37.6 agree with [17] 



-?'(r) = 



457r(2100027e 8 + 30148389005) 



137438953472 
confirming both (6) and our numerical accuracy. 



37.45 , (6) 
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FIG. 4: Plots of -iAG L (-«7). (a) For < 7 < 5.6. Insets 
show some regions with two close zeros, (b) The region 7 ~ 4 
expanded. The solid line is the numerical result, and the 
broken line is a fit to (8) with ai = —0.0227. 



The cut in G itself is 



AG(x,y; -ry) = -2179(7) 



f(x,-ij)f(y,-ij) 

J+(-h)J-(-h) 



(7) 



[cf. below (3) for the ± subscripts]. Although / can be 
stably integrated when matching to a Born approxima- 
tion instead of to /(£—>— 00, —17) ~ e~ 7Z [22], we prefer 
Jaffe's series [24, 27]. Figure 3 shows a typical result. For 
some refinements, see Ref. [23]. For 7 | 0, (7) reproduces 
G(s,y;t-K») [18]. 

The important limit is x, —y — ► 00. Since f(y, —i"f) ~ 
e 7 ^ and in general also / [x, — 17) ~ const • e lx , (7) 
has a strong position dependence — simply a result of 
the long propagation time, and removed if t is mea- 
sured from the first arrival at to(x,y) = x — y. Define 
G L (t') = lim^_^ 00 G(a;,y;to+t / ), so G L = J" 1 by (4) 
and the normalization of /, g. Thus, AG L = JI 1 — JZ 1 ; 
see Fig. 4 for results. 

The cut AG L (— ry) vanishes at (a) the zeros of 9(7) 
[cf. (7)], and (b) the anomalous points 7 = 5,1,5,..., 
where / and hence J diverge. The former depend only on 
V(x^>+oo) [cf. below (5)], the latter only on 00), 
scaling with A. If the two tails are separately adjusted 
the sequences are independent, but for A = 1 they share 



one member at T. Some members of the two sequences 
are close, cf. the insets in Fig. 4a. The many zeros also 
keep AG L small generically. 

Unconventional QNM. Although AG L (-zr) = 0, 
surprisingly AG L is largest (close to a dipole) near 7 = T 
(Fig. 4b; the dominance is less pronounced for finite y, cf. 
Fig. 3). This reveals a pair of nearby QNM poles u>±: if 
G±((j) w (ai±ia,2)/(iu)— 4— b±ic), with ai,a,2,b real and 
c > so uj± are on the unphysical sheets, then 



AG L H 7 ) 



2»a 2 ( 7 -r) 

( 7 -r-fe) 2 + c 2 



Here, a\C + a 2 b = enforces the zero at 7 = T. 
broken line in Fig. 4b shows this fit, yielding 

oj± + iT = =fc - ib w tO.027 + 0.0033i . 



(8) 



The 



(9) 



The extrapolation into the unphysical sheet can also 
be carried out analytically, by assuming that J + (w«r2) 
can be linearized up to the nearest zero, 



(10) 



Since AJ(O) = 0, the sheet of J(Q) need not be indi- 
cated. Following the methods of [17], one readily obtains 
J ( n ) = _ 5T?iul- Further, a lengthy calculation gives [23] 



49 2 38 h 17122265640585( ^E+ In8-i7r) 

- 245810518235861775 Ei(8+^) e - 8 
+ 36326230655979688] , (11) 



where Ei(z) = f* dt e l /t and r\ > is an infinitesimal. 
Insertion into (10) yields 



lu+ + iT w -0.03248 + 0.003436i ; 



(12) 



especially the agreement of Imw+ with (9) is remarkable. 
Since the latter value is also found by extrapolation, it 
need not be more accurate. 

Discussion. Like the hydrogen spectrum, the 
Schwarzschild spectrum contains discrete and contin- 
uum parts. While the former is a classic of physics, the 
latter is much more difficult, being spheroidal rather 
than hypergeometric [24]. We have characterized the 
continuum, recovering the behavior both for 7 J. and 
near the miraculous point T. This leads to the identifi- 
cation of an essentially new type of damped excitation 
lu± , which clearly affects the dynamics more than QNMs 
on the physical sheet at larger |Imw| [6]. An obvious 
question is whether there are further singularities on the 
unphysical sheets, and what is their influence on the cut. 

In a wider context, consider the Kerr hole. By com- 
paring numerics for moderately small a [16] with the 
QNM multiplct found analytically to branch off from fi 
at a = 0, one concludes that one additional multiplct has 
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to emerge (as a increases) near f2. Rather than the possi- 
bilities considered in Ref. [17], this multiplet may well be 
due to uj± splitting (since spherical symmetry is broken) 
and moving through the NIA as a is tuned. As a first 
step, the continuum should be studied also for a > 0. 

One is led to consider another Fourier contour go- 
ing into the unphysical sheet and detouring around u)± 
(Fig. 1, line b), including them as QNMs. This 
slightly reduces the continuum (often neglected as "back- 
ground") [23]. More importantly, if these poles emerge 
onto the physical sheet when tuning a parameter (say, a) , 
the total QNM and continuum contributions now become 
separately continuous. 

All of the above refers to £ = 2. The much larger 7 w 
T(£>3) are still unattainable numerically, but since any 
unconventional poles can be shown to be further away 
from [23], their influence on the cut should be smaller. 

These questions may be explored through solvable 
models with potential tails. Some aspects of the 
Regge- Wheeler equation can also be analyzed asymptot- 
ically [26]. Numerical algorithms valid on the NIA and 
even into the unphysical region would also be useful, al- 
lowing QNMs there to be studied directly rather than 
through extrapolation. 
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